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Abstract 

We examine the effects of the R parity odd renormalizable interactions on flavor changing rates 
and CP asymmetries in the production of fermion-antifermion pairs at leptonic (electron and muon) 
colliders. In the reactions, l~ fj+fji, [I — e, fi; J ^ J'] the produced fermions may be leptons, 
down-quarks or up-quarks, and the center of mass energies may range from the Z-boson pole up to 1000 
GeV. Off the Z-boson pole, the flavor changing rates are controlled by tree level amplitudes and the 
CP asymmetries by interference terms between tree and loop level amplitudes. At the Z-boson pole, 
both observables involve loop amplitudes. The lepton number violating interactions, associated with 
the couphng constants, Xijk, Kjkt ''■rs only taken into account. The consideration of loop amplitudes 
is restricted to the photon and Z-boson vertex corrections. We briefly review flavor violation physics 
at colliders. We present numerical results using a single, species and family independent, mass 
parameter, m, for all the scalar superpartners and considering simple assumptions for the family 
dependence of the R parity odd coupling constants. Finite non diagonal rates (CP asymmetries) entail 
non vanishing products of two (four) different coupling constants in different family configurations. 
For lepton pair production, the Z-boson decays branching ratios, Bjji = B{Z + Zj",), scale in 

order of magnitude as, Bjji ~ ((n")*(^^^^^^)'^'^ 10~^, with coupling constants A — Xijk or X[jk in 
appropriate family configurations. The corresponding results for d- and u-quarks are larger, due to 
an extra color factor, Nc — 3. The fiavor non diagonal rates, at energies well above the Z-boson pole, 
slowly decrease with the center of mass energy and scale with the mass parameter approximately as, 
ajji ~ ( oT )^ ~ ^(^ ~ W)f barns. Including the contributions from an sneutrino s-channel 
exchange could raise the rates for leptons or d-quarks by one order of magnitude. The CP-odd 
asymmetries at the Z-boson pole, Ajj' — b''''',^b''','' ? vary inside the range, (1Q~^ ^ 10^'^) sin 
where tp is the CP-odd phase. At energies higher than the Z-boson pole, CP-odd asymmetries for 
leptons, d-quarks and u-quarks pair production lie approximately at, (1Q~^ — 10~^) sin-i^, irrespective 



of whether one deals with light or heavy fiavors. 
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1 Introduction 



An approximate R parity symmetry could greatly enhance our insight into the supersymmetric flavor 
problem. As is known, the dimension four R parity odd superpotcntial trilinear in the quarks and leptons 
superfields, 

adds new dimensionless couplings in the family spaces of the quarks and leptons and their superpartners. 
Comparing with the analogous situation for the Higgs-meson-matter Yukawa interactions, one naturally 
expects the set of 45 dimensionless coupling constants, Xijk = —Xjik, Kjk^ Kjk ~ ~Kkj^ exhibit 
a non-trivial hierarchical structure in the families spaces. Our goal in this work will be to examine a 
particular class of tests at high energy colliders by which one could access a direct information on the 
family structure of these coupling constants. 

The R parity symmetry has inspired a vast literature since the pioneering period of the early 80's 
[§ |, |, |, |, I, 0, I and the maturation period of the late 80's and early 90's §, |l^, [l], 
This subject is currently witnessing a renewed interest |l6|. As is well known, the R parity odd 
interactions can contribute at tree level, by exchange of the scalar superpartners, to processes which 
violate the baryon and lepton numbers as well as the leptons and quarks flavors. The major part of the 
existing experimental constraints on coupling constants is formed from the indirect bounds gathered from 
the low energy phenomenology. Most often, these have been derived on the basis of the so-called single 
coupling hypothesis, where a single one of the coupling constants is assumed to dominate over all the 
others, so that each of the coupling constants contributes once at a time. Apart from a few isolated cases, 
the typical bounds derived under this assumption, assuming a linear dependence on the superpartner 
masses, are of order, [A, A', A"] < (10~^ - 10~^) looSey 

One important variant of the single coupling hypothesis can be defined by assuming that the dom- 
inance of single operators applies at the level of the gauge (current) basis fields rather than the mass 
eigenstate fields, as was implicit in the above original version. This appears as a more natural assumption 
in models where the presumed hierarchies in coupling constants originate from physics at higher scales 
(gauge, flavor, or strings). Flavor changing contributions may then be induced even when a single R 
parity odd coupling constant is assumed to dominate . While the redefined mass basis superpotcntial 
may then depend on the various unitary transformation matrices, V^'^, two distinguished predictive 
choices are those where the generation mixing is represented solely in terms of the CKM (Cabibbo- 
Kobayashi-Maskawa) matrix, with flavor changing effects appearing in either up-quarks or down-quarks 
flavors [r^ . A similar situation holds for leptons with respect to the couplings, Xijk, and transformations, 

A large set of constraints has also been obtained by applying an extended hypothesis of dominance of 
coupling constants by pairs (or more) . Several analyses dealing with hadron flavor changing effects (mixing 
parameters for the neutral light and heavy flavored mesons, rare mesons decays such as, K ^ it v + 
...) 0; lepton flavor changing effects (leptons decays, ^ + C + ' @ ^ conversion 

processes, neutrinos Majorana mass ...); lepton number violating effects (neutrinoless double 
beta decay [p2| |2^ , p^); or baryon number violating effects (proton decay partial branchings rare 
non-leptonic decays of heavy mesons nuclei desintegration |27|,...) have led to strong bounds on 
a large number of quadratic products of the coupling constants. All of the above low energy works, 
however, suffer from one or other form of model dependence, whether they rely on the consideration of 
loop diagrams |p5| , on additional assumptions concerning the flavor mixing |p^ , 19, or on hadronic 
wave functions inputs (26, 27 



Proceeding further with a linkage of R parity with physics beyond the standard model, our main 
observation in this work is that the R parity odd coupling constants could by themselves be an independent 
source of CP violation. Of course, the idea that the RPV interactions could act as a source of superweak 
CP violation is not a new one in the supersymmetry literature. The principal motivation is that, whether 
the RPV interactions operate by themselves or in association with the gauge interactions, by exploiting 
the absence of strong constraints on violations with respect to the flavors of quarks, leptons and the scalar 
superpartners by the RPV interactions, one could greatly enhance the potential for observability of CP 
violation. To our knowledge, one of the earliest discussion of this possibility is contained in ref.||^, where 
the role of a relative complex phase in a pair of A^^j, coupling constants was analyzed in connection with 
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the neutral K, K mesons mixing and decays and also with the neutron electric dipole moment. This 
subject has attracted increased interest in the recent literature |28[ |2^, ^ |3[ |3^, |3^, 

Thus, the role of complex A^^ j, coupling constants was considered in an analysis of the muon polarization 
in the decay, — > /i+ + + 7 [Q, and also of the neutral B, B meson CP-odd decays asymmetries 
p9| , |l[ ^ ; that of complex Xijk interactions was considered in a study of the spin-dependent asymmetries 
of sneutrino-antisneutrino resonant production of r— lepton pairs, ^ v^v ^ t^t~ and that of 
complex A^jj. interactions was considered as a possible explanation for the cosmological baryon asymmetry 
Q, as well as in the neutral B, B decays asymmetries An interesting alternative proposal @ 

is to embed the CP-odd phase in the scalar superpartner interactions corresponding to interactions of 
■^ijkKjk type. Furthermore, even if one assumes that the R parity odd interactions are CP conserving, 
these could still lead, in combination with the other possible sources of complex phases in the minimal 
supersymmetric standard model, to new tests of CP violation. Thus, in the hypothesis of pair of dominant 
coupling constants new contributions involving the coupling constants A^^^, and the CKM complex phase 
can arise for CP-odd observables associated with the neutral mesons mixing parameters and decays [ p9| , 
^ Also, through the interference with the extra CP-odd phases present in the soft supersymmetry 
parameters. A, the interactions Xijk and A^^^ may induce new contributions to electric dipole moments 

We propose in this work to examine the effect that R parity odd CP violating interactions could 
have on flavor non-diagonal rates and CP asymmetries in the production at high energy colliders of 
fermion-antifermion pairs of different families. We consider the two-body reactions, l~{k) + l^{k') 
fj{p) + Ij'{p')i [J 7^ J'] where I stands for electron or muon, the produced fermions are leptons, down- 
quarks or up-quarks and the center of mass energies span the relevant range of existing and planned 
leptonic (electron or muon) colliders, namely, from the Z-boson pole up to 1000 GeV. High energy colliders 
tests of the RPV contributions to the flavor diagonal reactions were recently examined in ^ 
and for flavor non-diagonal reactions in . 

The physics of CP non conservation at high energy colliders has motivated a wide variety of proposals 
in the past |Q and is currently the focus of important activity. In this work we shall limit ourselves to the 
simplest kind of observable, namely, the spin independent observable involving differences in rates between 
a given flavor non-diagonal process and its CP conjugated process. While the R parity odd interactions 
contribute to flavor changing amplitudes already at tree level, their contribution to spin independent CP- 
odd observables entails the consideration of loop diagrams. Thus, the CP asymmetries in the Z-boson 
pole branching fractions, B{Z —^fj + fj'), are controlled by a complex phase interference between non- 
diagonal flavor contributions to loop amplitudes, whereas the off Z-boson pole asymmetries are controlled 
instead by a complex phase interference between tree and loop amplitudes. Finite contributions at tree 
level order can arise for spin dependent CP-odd observables, as discussed in refs. |3^ . 

It is useful to recall at this point that contributions in the standard model to the flavor changing 
rates and/or CP asymmetries can only appear through loop diagrams involving the quarks-gauge bosons 
interactions. Corresponding contributions involving squarks-gauginos or sleptons-gauginos interactions 
also arise in the minimal supersymmetric standard model. In studies performed some time ago within the 
standard model, the flavor non diagonal vector bosons (Z-boson and/or W-bosons) decay rates asymme- 
tries H H, and CP-odd asymmetries MM were found to be exceedingly small. (Similar conclusions 
were reached in top-quark phenomenology ]5C|.) On the other hand, in most proposals of physics beyond 
the standard model, the prospects for observing flavor changing effects in rates ^ 48, ^ 

or in CP asymmetries Q, |5^ are on the optimistic side. Large effects were reported for the supersym- 
metric corrections in flavor changing Z-boson decay rates arising from squarks flavor mixings [p3[|, but 
the conclusions from this initial work have been challenged in a subsequent work [ p4| involving a more 
complete calculation. 

The possibility that the R parity odd interactions could contribute to the CP asymmetries at observ- 
able levels depends in the first place on the accompanying mechanisms responsible for the flavor changing 
rates. Our working assumption in this work will be that the R parity odd interactions are the dominant 
contributors to flavor non-diagonal amplitudes. 

The contents of this paper are organized into 4 sections. In Section |^, we develop the basic formalism 
for describing the scattering amplitudes at tree and one- loop leve ls. We discuss the case of leptons, 
down-quarks and up-quarks successively in subsection s 2.2 and 2.3. The evaluation of the one-loop 



loop diagrams is based on the standard formalism of |£6|. Our calculations here closely parallel similar 
ones developed ||5^, ^ in connection with corrections to the Z-boson partial widths. In Section ^, we 
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first briefly review the physics of flavor violation and next present our numerical results for the integrated 
cross sections (rates) and CP asymmetries for fermion pair production at and off the Z-boson pole. In 
Section ^, we state our main conclusions and discuss the impact of our results on possible experimental 
measurements. 

2 Production of fermion pairs of different flavors 

In this section we shall examine the contributions induced by the RPV (R parity violating) couplings 
on the flavor changing processes, l^{k) + l^{k') fj{p) + fj'{p'), [I — e, fi; J ^ J'] where / stands 
for leptons or quarks and J, J' are family indices. The relevant tree and one- loop level diagrams are 
shown schematically in Fig. ^ At one-loop order, there arise 7— and Z— boson exchange triangle 



i(k) fj(P) 




(a) (b) (c) 



Figure 1: Flavor non-diagonal process of l^l^ production of a fermion- antifcrmion pair, l^{k) + l^{k') 
fj{p) + fj'{p')- The tree level diagrams in (a) represent t- and s-channel exchange amplitudes. The loop 
level diagrams represent 7 and Z gauge boson exchange amplitudes with dressed vertices in (b) and box 
amplitudes in (c). 

diagrams as well as box diagrams. In the sequel, for clarity, we shall present the formalism for the 
one-loop contributions only for the dressed Zff vertex in the Z-boson exchange amplitude. The dressed 
7-exchange amplitude has a similar structure and will be added in together with the Z-boson exchange 
amplitude at the level of the numerical results. Since we shall repeatedly refer in the text to the R 
parity odd effective Lagrangian for the fermions-sfermion Yukawa interactions, we quote below its full 
expression, 

i = ^ j ^Ay-fe [C'tLekRBjL + ejLekRViL + e*kR^iRf^jL ])] 
ijk ^ _ 

+ Kjkl^iLdkRdjL + djLdkRViL + dXjij^iRdjL — CiLdkRUjL — UjLdkRBiL — dXLlC^j^UjL] 

+ i3k<^aP^[u*^Rdof}Rdl^L + d*jpRUtaRdl^L + dt^R^zaRd^jpL " U ^ ^)]| + ^.C. , (2) 

noting that the summations run over the (quarks and leptons) families indices, i, j, k = [(e, /i, t); (d, s, b); (u, c, t)], 
subject to the antisymmetry properties, A^-fe = —\jik, ^'Ijk = ~^ikj- precedence conventions for 

operations on Dirac spinors such that charge conjugation acts first, chirality projection second and Dirac 
bar third, so that, i'l r = {'^'^)l.r- 

2.1 Charged lepton-antilepton pairs 
2.1.1 General formalism 

The process l~{k) + l^{k') &j{p) + ej'(p')j for I = e, n; J ^ J' , can pick up a finite contribution at 
tree level from the R parity odd couplings, Xijk, only. For clarity, we treat in the following the case of 
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electron colliders, noting that the case of muon colliders is easily deduced by replacing all occurrences in 
the RPV coupling constants of the index 1 by the index 2. There occur both t-channel and s-channel 
ViL exchange contributions, of the type shown by the Fcynman diagrams in (a) of Fig. |^. The scattering 
amplitude at tree level, Mt, reads: 
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' \ll\uj,Vii{k')uL{k)uL{p)vii{p') + \*ii\iJ,J'VL{k')uB.{k)uB{p)vL{p') 



(3) 



where to obtain the saturation structure in the Dirac spinors indices for the t-channel terms, we have 
applied the Fierz rearrangement formula, uji{p)uL{k)v L{k')v ii{p') = \uR{p)^^j,VR[p')vL{k')^^UL{k). The 
t-channel (s-channel) exchange terms on the right hand side of eq.(^) include two terms each, called TZ- 
and £-type, respectively. These two terms differ by a chirality flip, L R, and correspond to the distinct 
diagrams where the exchanged sneutrino is emitted or absorbed at the upper (right-handed) vertex. 
The Z-boson exchange amplitude (diagram {b) in Fig. |l|) at loop level. Mi, reads: 

2 



M, 



J J' 



2 cos 9 



w 



v{k')-iAa[eL)PL+a{eii)PR ]u{k) 



1 



(4) 



where the Z-boson current amplitude vertex function, T^{p,p'), is defined through the effective La- 
grangian density. 



L 



2 cos bw 

For later convenience, we record for the processes, Z{P = p + p') f{p) + f'{p') and Z{P) fnip) + 
f^{p'), the familiar definitions of the Z-boson bare vertex functions. 



/(p)7m aifL)PL + a{fn)Pn /'(p') + {p - p')MP>UH)fH{p) 



(5) 



where the quantities denoted, a^fn) = an if) and aifn), taking equal values for both fermions and 

2Qxw, where H = (L, R), xw = sin^ 9w, are 



sfermions, are defined by, aifu) = a(/ff) = "^T^if) 
SU {2)h Cartan subalgebra generators, and Q = + Y, Y are electric charge and weak hypercharge. 
These parameters satisfy the useful relations: a[f^) = —a{fH), ol(/'^) = —aii{f), ajiif^) = —aL{f). 

Throughout this paper we shall use the conventions in Haber-Kane review |]59| (metric signature (H ), 

P^t-j = (1 =p 75)72, etc..) and adopt the familiar summation convention on dummy indices. 

The Lorentz covariant structure of the dressed Z-boson current amplitude in the process, Z{P) —f 
fj{p) + fj'{p')i for a generic value of the Z-boson invariant mass s = P^, involves three pairs of vectorial 
and tensorial vertex functions, which are defined in terms of the general decomposition: 



Tlip.p') = u{p) 



7m 



A 



-a 



mj + mj' 



t'\f)PL 
t^u\{p + 



A-'/{f)PR 



+ 75d'"'] + (p-p')''[*^''''' +756 



(6) 



where, cr^^ = §[7^,7;^]. The vector vertex functions separate additively into the classical (bare) and 
loop contributions, Aj/ (/) = aH{f)Sjj' + Aj/ (/), [H = L,R]. The tensor vertex functions, associated 
with a'^'^{p + p')i/, include the familiar magnetic and electric Z f f couplings, such that the flavor 

:[a'^'^, d'^'^], identify, in the small momentum transfer limit, with 



diagonal vertex functions, —75 — ^-3— „^ 

o ' 2 cos &w 2m ,7 

the fermions Z-boson current magnetic and (P and CP-odd) electric dipole moments, respectively. In 
working with the spinors matrix elements, it is helpful to recall the mass shell relations, u{p)^ 
= mju{p), '^'v{p') — ~mjiv{p'), and the Gordon type identities, appropriate to the saturation of the 
Dirac spinor indices, u{p) ■ ■ ■ v{p'), 



It 



ia^uipTP . 



(mj + TO,/')7^ 



75 
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Based on these identities, one also checks that the additional vertex functions, [b'^"^ , e'^"^ ], associated 
with the Lorentz covariants, a^^{p — 75], can be expressed as linear combinations of the vector or 

axial covariants, 7^ [1,75], and the total momentum covariants, (p + p')^ [1j75]- The latter will yield, 
upon contraction with the initial state Zl~l~^ vertex function, to negligible mass terms in the initial 
leptons. 

Let us now perform the summation over the initial and final states polarizations for the summed 
tree and loop amplitudes, Af = M/^' + Mf"^ , where the lower suffices t, I stand for tree and loop, 
respectively. (We shall not be interested in this work in spin observables.) A straightforward calculation, 
carried out for the squared sum of the tree and loop amplitudes, yields the result (a useful textbook to 
consult here is ref. Q): 



E 

pol 



M, 



JJ'\2 
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'a{eL)A-j/ (e,s + ie) 



2{t-ml^J ' V2cose'H' 
+ 8mjmj,(k ■ k')^LL{n)+8ml{p-p')ipnn{TZ) + 

'^a{en)Al'''{e,s + ie) ^ 



rrir, 



imzT2 



lQ{k ■ p)(k' ■ p') 



9_ 

2 cos e 



w 



imzT2 



l&{k ■ p){k' ■ p') 



+ 8mjmj'{k ■ k')(pRR{C) + 8mg(p • p')ipLL{C) + 8 



{k ■ k'){p-p), 



(7) 



where we have introduced the following functions, associated with the TZ- and £-type contributions: 



Vhh' (>C) = 



g 



a{eH)A-l/,' {e,s + ie)Y f ^iuKu' 



2 cos 6w J V s — m'^ + imzTz J \2{t~mi,^) 
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2 cos 9w 



imzTz 



(8) 



The two sets of terms in eqs.(^) and labelled by the letters, TZ, C, are associated with the two 
t-channel exchange contributions in the tree amplitude, eq.(|^), which differ by the spinors chirality 
structure and the substitutions, AujA^^j, — > X*j^Xijri. The terminology, £, TZ, is motivated by the fact 
that these contributions are controlled by the Z-boson left and right chirality vertex functions, A^ and 
Afj, respectively, in the massless limit. 

The imaginary shift in the argument, s + ie (representing the upper lip of the cut real axis in the 
complex s-plane) of the vertex functions, Aj/ (/, s + ie), has been appended to remind us that the one- 
loop vertex functions are complex functions in the complex plane of the virtual Z-boson mass squared, 
s = {p+p'Y' , with branch cuts starting at the physical thresholds where the production processes, such as, 
Z^f + foiZ^f + f*, are raised on-shell. For notational simplicity, we have omitted writing several 
terms proportional to the initial leptons masses and also some of the small subleading terms arising from 
the loop amplitude squared. At the energies of interest, whose scale is set by the initial center of mass 
energy or by the Z-boson mass, the terms involving factors of the initial leptons masses nie, are entirely 
negligible, of course. Thus, the contributions associated with iprr{TZ), lpll{L) can safely be dropped. 
Also, the contribution from ifLLiTZ) and i^rr^C) which are proportional to the final state leptons masses, 
mj, and mji, can to a good approximation be neglected for leptons production. Always in the same 
approximation, we find also that interference terms are absent between the s-channel exchange and the 
t-channel amplitudes and between the s-channel tree and Z-boson exchange loop amplitudes. Similarly, 
because of the opposite chirality structure of the first two terms in M/'^ , their cross-product contributions 
give negligibly small mass terms. 



2.1.2 CP asymmetries 

Our main concern in this work bears on the comparison of the pair of CP conjugate reactions, l^{k) + 
l^{k') ey(p) + e^j,{p') and l~ {k) + l'^{k') e^r{p) + £^{p')- Denoting the summed tree and one-loop 
probabihty amplitudes for these reactions as, M'^-'' = M^-^' + Mf-^' , M'^-'' = M/'"' -I- Mf •' = M-''-^ , 
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we observe that these amphtudes are simply related to one another by means of a specific complex 
conjugation operation. The general structure of this relationship can be expressed schematically as: 



M 



J J' 



"0 



-,JJ' ipJJ 



J J' 



-,.7,7'* 



(9) 



where for each of the equations above, referring to amplitudes for pairs of CP conjugate processes, 
the first and second terms correspond to the tree and loop level contributions, with Cq"^ — ^.JJ * 

representing the tree amplitudes and a^'^ , a;^ ^ 



andF, 



J J' 



F, 



J' J 



F^'^ representing the complex 
valued coupling constants products and momentum integrals in the loop amplitudes. The functions 
pjj Ysmgi symmetric under the interchange, J <-!■ J'. The summation index a labels the family 
configurations for the intermediate fermions-sfermions which can run inside the loops. Defining the CP 
asymmetries by the normalized differences. 



\M 



J J' 12 



\M 



JJ'\2 



and inserting the decompositions in cq. 



Ajj' 



|ao| 



\M-'J'\^ + \MJJ'\2' 
, the result separates additively into two types of terms: 

^ Im{aoa*^)Im{Fa{s + ie)) 



Im{aaa^,)Im{Fa{s + ie)F*,{s + ie)) 



a<a' 



(10) 



J J' „JJ' 



and 



where, for notational simplicity, we have suppressed the fixed external family indices on a, 
F^^ , and replaced the full denominator by the tree level amplitude, since this is expected to dominate 
over the loop amplitude. The first term in ( p^ is associated with an interference between tree and loop 
amplitudes and the second with an interference between terms arising from different family contributions 
in the loop amplitude. In the second term of eq.(|l^), the two imaginary parts factors are antisymmetric 
under the interchange of indices, a and a', so that their product is symmetric and allows one to write, 
^a<a' ^ 5 X^a^a' • ^o obtain a more explicit formula, let us specialize to the specific case where 
the Z-boson vertex functions decompose as, Aj/ (/, s + ie) = J2a^^J'^Ha + *^)- The first factors, 
^JJ' — ^ijjKjj' (using a — (ij) and notations for the one- loop contributions to be described in the next 
subsection), include the CP-odd phase from the R parity odd coupling constants. The second factors, 
^Ha ' include the CP-even phase from the unitarity cuts associated to the physical on-shell intermediate 
states. In the notations of eq.(0). 
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)'a{eH')b1y,, Fi'' - FIj'^ {s + ie)/{s - m% + imzTz), 



-J7' 



2 cos dw 



where the right hand sides incorporate appropriate sums over the chirality indices, H' , H oi the initial 
and final fermions, respectively. 

Applying eq.(^) to the asymmetry integrated with respect to the scattering angle, one derives for the 
corresponding integrated tree-loop interference contribution. 



A 



J7' 



= -4 



9 



2 cos 6w 

dx tttt: 2 



(1 - x)^ 



J2\^^lJKlJ'\' f dx 
J — 1 



(1-x)^ 



4(i - )2 



imzT2 



(11) 



where, 9, [x — cos 9] denotes the scattering angle variable in the center of mass frame and the Mandelstam 



variables in the case of massless final state fermions take the simplified expressions, s = {k + k') , t 



2 + - 



{k-pY 



3(1 — X), U : 



final fermions masses, mj, 



-={k- 



si- 



-P) - 2 
are: -^/s = 2k 



Useful kinematical relations in the general case with 



where, E„ 



is 



Ep + Epi, t 
2 



'sEp{l-l3x) 



m2,)/(2yi), Ep, = {s + m^j, - m^j)/{2^s), /3 = p/Ep, /?' 



= P/Ep. 



k, p denoting the center of mass momenta of the two-body initial and final states, respectively, 
unpolarized differential cross section reads then, da/dx = i28ns\k\ Spo/ l-^^P- 



-P'x), 
with 
The 
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For the Z-boson pole observables, the flavor non-diagonal branching ratios and CP asymmetries (where 
one sets, s — m|) are defined in terms of the notations specified in the preceeding paragraph by the 
equations, 

R _ T{Z^fj + fj,)+T{Z^f,r+fj) JAiJ'if)? + \Al,''if)\' 



A 



J J' 



nz ^fj+ ir) - nz ^ fj- + fj) 

T{Z ^ fj + fj')+ T{Z fj, + fj) 

J:H=L,R\J:.bf^,{f)Fg{s + ^e)\^ " ^ ^ 

For completeness, we recall the formula for the Z-boson decay width in fermion pairs (masslcss limit), 

r{z ^ fj + 1,,) = ^I^[\Ai^\f)\^ + 

where, c/ = [1, Nc], for [/ — l,q] {Nc = 3 is the number of colors in the SU{3)c color group) and the 
experimental value for the total width. T{Z — > all)exp = 2.497 GeV. 

The expressions in eqs.(pl|) and ( [l^ ) for the CP asymmetries explicitly incorporate the property of 
these observables of depending on combinations of the RPV coupling constants, such as, Arg{XiijX*^j, X*,jjXiijj>), 
or Arg{Xij J \*j J, X*, J, J Xi> ji ji) , which are invariant under complex phase redefinitions of the fields. This 
freedom under rephasings of the quarks and leptons superfields actually removes 21 complex phases from 
the complete general set of 45 complex RPV coupling constants. 

2.1.3 One-loop amplitudes 

The relevant triangle Feynman diagrams, which contribute to the dressed Z-boson leptonic vertex, 
Z(P)/^(p)Z+(p'), appear in three types, fermionic, scalar and self-energy, as shown in Fig. |[ We 
consider first the contributions induced by the R parity odd couplings, A^^^. The intermediate lines can 




Figure 2: One-loop diagrams for the dressed Z{P) f{p)f{p') vertex. The flow of four-momenta for the 
intermediate fermions in (a) is denoted as, Z{P) f{Q) + f{Q') fj{p) + fj'{p')- Similar notations are 
used for the sfermions diagram in (b) where, Z{P) f'{Q) + f *{Q'), and for the self-energy diagrams 
in (c). 

assume two distinct configurations which contribute both, in the limit of vanishing external fermions 
masses, to the left-chirality vertex functions only. We shall refer to such contributions by the name C- 
type contributions, reserving the name 7?.- type to contributions to the right-chirality vector couplings. 
The two allowed configurations for the internal fermions and sfermions are: / = ('^c);f' = (~'^)- Our 

calculations of the triangle diagrams employ the kinematical conventions for the flow of electric charge 
and momenta indicated in Fig.0, where P = p + p' = Q + Q' — k + k' . The summed fermion and scalar 
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Z-boson current contributions are given by: 



(g^-mp((Q-p- 



P. 



m),){{Q-pY 



VOj) 



(13) 



The integration measure is defined as, Jg = (2^^ / For ^ convenient derivation of tlie self-energy 
diagrams, one may invoke the on shell renormalization condition which relates these to the fields renor- 
malization constants. Defining schematically the self-energy vertex functions for a Dirac fermion field ip 
by the Lagrangian density, L = i'tp{^— m + 'E{p))ip, = mao + '^{a'^PL + o'^Pr), the transition from 

bare to renormalized fields and mass terms may be efltected by the replacements, 



[H = L, R] 



By a straightforward generalization to the case of several fields, labelled by a family index J, the fields 
renormalization constants become matrices, Zfj, = (1 + a^)^j,. The self-energy contributions to the 
dressed Z-boson vertex function is then described as, 



^Up,p')se 



(iZ?j'Z'/,y)i l)uiphMfH)PHvip') 

H=L,R ^ ^ 

- E \{'^lr{p) + ^'}',{p'mp)l,a{fH)PHv{p'), 



H=L,R 



where for the case at hand, 

S,7J'(P) 



mf)PL 



(-Q2 + m2)(-(g-p)2 + m|,)^ 



(14) 



(15) 



so that Ojj, = and (Tq = 0. Similar Feynman graphs to those of Fig. ||, and similar formulas to those 
of eqs.(13) and (p^, obtain for the dressed photon current case, 7(P)Z~(p)Z+(p'). 

We organize our one- loop calculations in line with the approach developed by 't Hooft and Veltman [ p5[ 
and Passarino and Veltman [Q, keeping in mind that our spacetime metric has an opposite signature to 
theirs, ( — \--\ — h). For definiteness, we recall the conventional notations for the two-point and three-point 
integrals. 



(2^) 



-[Bo,-p^Si] - / 

JQ 



{~Q^+ml){-{Q-pr + miy 



(16) 



[Co, -P,j.Cll-Pf,Ci2, Pf,PuC21+plp'f,C22 + {Pf,P^+PuPf,)C23 - 9,^1,024] 

{-Q^ + ml){-{Q-py + ml)i-iQ-p-p')2 + ml)' ^ ' 

where the arguments for the B— and C— functions are defined as: Ba{—p, mi, 7712), [A — 0,1] and 
Cb{~P, —p' , fni, 1712, W3), [B — 0, 11, 12, 21, 22, 23, 24]. In the algebraic derivation of the one-loop am- 
plitudes, we find it convenient to introduce the definitions: = + Pf^, p'^ = — p^, where 
P = p + p', p = ^{p — p'). The terms proportional to the Lorentz covariant — {p + p')^ will then 
reduce, for the full Z-boson exchange amplitude in eq. (^), to negligible mass terms in the initial leptons. 

Dropping mass terms for all external fermions, the tensorial couplings cancel out and we need keep 
track of the vector couplings only, with the result: 



aihWfCo + a{fR) (^4^) - 2C24 - mj,Co^ + 2a(/i)C24 + ^(/l)^^ 



A'r'\C) = 0. (18) 
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The cancellation of the right chirality vertex function in this case is the reason behind our naming 
these contributions as £-type. The two-point and three-point integrals functions without a tilde symbol 
arise through the fermion current triangle contribution and the self-energy contribution (represented 

(2) 

by the term proportional to B\ ). These involve the argument variables according to the following 
conventions, b'^^^ ~ Ba{^P — p' ,rnf,mf), B^^ = BA{—p,iTif,mp), B^^^ = Ba{~p' ,rn^,,mf), 

and Cb — CA{—p,—p'T''nf,mp,mf). The integral functions with a tilde arise in the sfermion current 

diagram and are described by the argument variables, Ca = Ca{—P, —p'^rnp, nif^mp). 

A very useful check on the above results concerns the cancellation of ultraviolet divergencies. This 
is indeed expected on the basis of the general rule that those interaction terms which are absent from 
the classical action, as is the case for the flavor changing currents, cannot undergo renormalization. A 
detailed discussion of this property is developed in [ |6l| . The logarithmically divergent terms in eq. ( p^ ) , 
proportional to the quantity, A = —j^zi + 7 — InTr, as defined in |56j, arise from the two- and three- 
point integrals as, Bq A, Si — > ^"24 |A, all other integrals being finite. Performing these 
substitutions, we indeed find that A comes accompanied by the overall factors, [— a(eL) -l-a(u^) + a{df()], 
or [—a(eL) + a^dfi) + aR{u'^)], which both do vanish in the relevant configurations for /, /'. 

Let us now consider the R parity odd Yukawa interactions involving the Xijk ■ These contribute through 
the same triangle diagrams as in Fig. ^. There arise contributions of £-type, in the single configuration, 
/ = Cfej/' = i^iL of 7?.- type in the two configurations, / = (JJ,^),/' = (g'^)- Following the same 
derivation as above, and neglecting all of the external mass terms, we obtain the following results for the 
one-loop vector coupling vertex functions: 

^J,7'(£) = ^^l^p HfLWfCo + a{fn)(B^o'^ - 2C24 - mj,Co) + 2a{fl)C2A + a(/L)i?P^], 
^''/(T^) = ^^J^^HhWfCo + a{fL){B'i^ - 2C24 - m),C^) + 2a(/i)C'24 + a(/7?)sf 

(19) 

with Aj^ (C) = 0, Aj/ (TZ) = 0. We note that the C, TZ contributions are related by a mere chirality 
flip transformation and that the color factor, TVc, is absent in the present case. 



2.2 Down-quark-antiquark pairs 

The processes involving flavor non-diagonal flnal down-quark-antiquark pairs, l~ {k) -f l^{k') dj{p) + 
dji{p'), pick up non vanishing contributions only from the A^^-^, interactions. Our discussion here will be 
brief since this case is formally similar to the leptonic case treated in subsection In particular, the 
external fermions masses, for all three families, can be neglected to a good approximation at the energy 
scales of interest. The tree level amplitude comprises an 7?.- type single t-channel w-squark exchange 
diagram and two s-channel diagrams involving v and v sneutrinos of the type shown in (a) of Fig. |^, 

'r—:UR{p)^^VR{p')vL{k')"fp.UL{k) 

KiiKj.j'VR{k')uL{k)uL{p)vR{p') + \*-^^\[jj,VL{k')uR{k)uR{p)vL{p') , 

(20) 

where a Kronecker symbol factor, Sab, expressing the dependence on the final state quarks color indices, 
d°'db, has been suppressed. This dependence will induce in the analog of the formula in eq.(^ expressing 
the rates, an extra color factor, N^- 

At one-loop level, the dressed Z dj dji vertex functions in the Z-boson s-channel exchange amplitude 
can be described by the same type of triangle diagrams as in Fig. ^. The fields configurations circulating 
in the loop correspond now to quarks-sleptons of £-type, / = d^; Z' = v*^, and of 7?.- type, / = (^^); I' = 

(g*^). There also occurs corresponding leptons-squarks fields configurations of £-type, / = ^f; /' = d^R, 

and 7?.— type, / — (^'); /' — (~^'^). The £— and TZ— type contributions differ by a chirality flip, the flrst 

contributing to A'['' and the second to Aj^^ . The calculations are formally similar to those in subsection 



JJ' 



2{t - m: 
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2.1 and the final results have a nearly identical structure to those given in (|T8|). For clarity, we quote the 
final formulas for the one- loop vector coupling vertex functions, 

^jj'(^) = ^'^ll^p HhWfC^ + a{fR){B^'^ - 2C24 - mj,Co) + 2a(/i)C24 + a(rfL)i3p'], 



(21) 



where the intermediate fermion-sfermion fields are labelled by the indices /, /'. There are implicit sums 
in eq.(|2l|) over the above quoted leptons-squarks and quarks-sleptons configurations. The attendant 
ultraviolet divergencies are accompanied again with vanishing factors, a{dii) — a{dL)~\-a{i''jn) = 0, a{dL) — 
a{dR) + a{vL) = 0. 



2.3 Up-quark-antiquark pairs 

The production processes of up-quark-antiquark pairs of different families, l~{k) + l~^{k') —> uj{p) + 
^J'{p')i niay be controlled by the A'^j, interactions only. The tree amplitude is associated with an u- 
channel J-squark exchange, of type similar to that shown by (a) in Fig. |l|, and can be expressed as. 



2{u 



) 



(22) 



after using the Fierz reordering identity, appropriate to commuting Dirac (rather than anticommuting 
field) spinors, u'^{k)PLv{p')u{p)Piiv'^{k') = +^VL{k')j^^UL{k)uL{p)j^VL{p')- We have omitted the Kro- 
necker symbol Sab on the u°'Ub color indices, which will result in an extra color factor Nc = 3 for the rates, 
as shown explicitly in eq.(^3|) below. The present case is formally similar to the leptonic case treated in 
subsection 2T, except for a chirality flip in the final fermions. We are especially interested here in final 
states containing a top-quark, such as tc or tu, for which external particles mass terms cannot obviously 
be ignored. The equation, analogous to (|^), which expresses the summations over the initial and final 
polarizations in the total (tree and loop) amplitude, takes now the form. 



E 

pol 



IJk 



'■a{eL)Ai'''{u,s + ie) 



2(u — m^- ) ' \2 cos 9w 
X 16(fc • p'){k' ■ p) + %mjmj. [k ■ k')^LR{C) 



imzTz 



where 0{m?^) terms were ignored and we have denoted. 



2 cos 9w 



imzTz 



V \ * 



2{u — 

^ dkR 



-\- C. C 



(23) 



(24) 



The modified structure for the kinematical factors in the above up-quarks case, eq.(23), in comparison 
with the leptons and d-quarks case, eq.(|^), reflects the difference in chiral structure for the RPV tree 
level amplitude. 

In the massless limit for both the initial and flnal fermions (where helicity, h = (—1, -1-1), and chirality, 
H = {L, R), coincide) the RPV interactions contribute to the helicity amplitudes for the process, — > 
/,/ + fj'i in the mixed type helicity configurations, hi- - 



(same as for the RFC 

gauge interactions) which are further restricted by the conditions, hi- = —hf,, for leptons and d-quarks 
production, and hi- = hfj, for up-quarks production. The dependence of the RFV scattering amplitudes 
on scattering angle has a kinematical factor in the numerator of the form, [l-\-hi-hfj cos6]^. [The parts in 
our formulas in eqs. ( p3| ) and (^, containing the interference terms between RFV and RFC contributions, 
partially agree with the published results 41 . We disagree with ^ on the relative signs of RFV 
and RFC contributions and with |^ on the helicity structure for the up-quarks case. Concerning the 
latter up-quarks case, our results concur with those reported in a recent study p3[.] 

The states in the internal loops of the triangle diagrams occur in two distinct £-type configurations. 



The calculations involved in keeping track of the mass terms are rather tedious. 
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They were performed by means of the mathematica software package, "Tracer" whose results were 
checked against those obtained by means of "FeynCalc" [Q. The relevant formulas for the vertex 
functions read: 



(47r)2 



aLiu) ) + a{fL)mf'Co + a{f) ( 2C24 + 2 (C12 - C21 + C23 - Cn ) 



( B'i^ - 2 C24 - m^.,2Co + m/ ( Co + SCn - 2Ci2 + 2C21 - 2C23 ) - m'j,Ci2 



(47r)2 



mjmji 



C23 + C; 



22 



Cii + C12 — 2C23 + 2C; 



22 



(25) 



The above formulas include an implicit sum over the two allowed configurations for the internal fermion- 
sfermions, namely, a{dkH), o,{e*]^) and a(e^^), aidkR). For completeness, we also display the formula 
expressing the tensorial covariants. 



mjPL aifR){Cu - C12 + C21 - C2 



fj. \Ft P Jtensorial — (47r)2 "^M"" 

a(/')(C'ii + C21 - (7i2 - C23)) + mj,Pn{ + a{fn){C22 - C23) + a(/')(C23 - C22) 



(26) 



The complete Zfjfj> vertex function, = T^g^t^riai + rJ;'e„sormi' should (after extracting the external 
Dirac spinors and the RPV coupling constant factors) be symmetrical under the interchange, J <-^- J', 
or more specifically, under the interchange, mj <-> mji. This property is not explicit on the expressions 
in eqs. (|2^) and (|2^), but can be established by reexpressing the Lorentz covariants by means of the 
Gordon identity. The naive use of eq.(p^) to compute CP-odd asymmetries would seem to yield finite 
contributions ^even in the absence of a CP-odd phase) from the mass terms in the vectorial vertex 
functions, Aj/ , owing to their lack of symmetry under, mj ^ mji. Clearly, this cannot hold true and 
is an artefact of restricting to the vectorial couplings. Including the tensorial couplings is necessary for 
a consistent treatment of the contributions depending on the external fermions masses. Nevertheless, we 
emphasize that the tensorial vertex contributions will not included in our numerical results. 

Finally, we add a general comment concerning the photon vertex functions, A^ ^j^ , and the way to 



incorporate the 7-exchange contributions in the total amplitudes, eqs.([7|) and (E3[). One needs to add 



terms obtained by substituting. 



e _ g sm Wyy 

2cosev^ ' 2 - — 2 — ' "^.wv./ ) ^vvj - "^z - 
along with the substitution of Z-boson by photon vertex functions, Ajfj^{e,s + ie) — > Aj^'j^ [e,s + le] 
The substitution which adds in both Z-boson and photon exchange contributions reads explicitly: 



aL,fl(/)^2Q(/), {s-ml + imzTzY 



^1% (e, s 



[aR,Li<^)Ai-^i^] 



aR,L{e) ^ a{f)Cf + 2Q{e) sin^ 0w cos^ Ow[{s - rn^z + 
/ 



Tz)ls\Y^2Q{f)Cf 
f 



where we have used the schematic representation, A'l'-'^^ — J2f "(/)C'/. 



3 Basic assumptions and results 

3.1 General context of flavor changing physics 

To place the discussion of the RPV effects in perspective, we briefly review the current situation of flavor 
changing physics. In the standard model, non-diagonal effects with respect to the quarks flavor arise 
through loop diagrams. The typical structure of one-loop contributions to, say, the Zff vertex function, 
Si ^Lj^iJ'H''^i^ /''^'z)^ involves a summation over quark famihes of CKM matrices factors times a loop 
integral. This schematic formula shows explicitly how the CKM matrix unitarity, along with the near 
quarks masses degeneracies relative to the Z-boson mass scale (valid for all quarks with the exception 
of the top-quark) strongly suppresses flavor changing effects. Indeed, for the down-quark-antiquark 
case, the Z-boson decays branching fractions, Bjj', were estimated at the values, 10~^ for (bs + sb), 
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10 ^ for {bd + db), 10 for (sd + ds), and the corresponding CP asymmetries, Ajj', at the values, 
[10-5 ,10-3, 10-i]sin(5cKM @, |4| , respectively. 

By contrast, flavor changing effects are expected to attain observable levels in several extensions of 
the standard model. Thus, one to three order of magnitudes can be gained on rates Bjji in models 
accommodating a fourth quark family |48| For the two Higgs doublets extended standard model, a 
recent comprehensive study of fermion-antifermion pair production at leptonic colliders |]5l| quotes for 
the flavor changing rates, Bjj, w IQ-'^ - 10"^ for Z ^ {b + s) + {s + b) and ajj, « 10"^ - lO^'^i?, where, 
R = cr(e+ + ^ /i+ + fi^) — 47ra^/(3s) — 86.8/ {y^y f barns {TeV)~^. Large CP violation signals are 
also found in the reaction, pp — tbX, in the two Higgs doublets and supersymmetric models 

For the minimal supersymmetric standard model, due to the expected nearness of superpartners 
masses to mz, flavor changing loop corrections can become threateningly large, unless their contribu- 
tions are bounded by postulating either a degeneracy of the soft supersymmetry breaking scalars masses 
parameters or an alignment of the fermion and scalar superpartners current-mass bases transformation 
matrices. An early calculation of the contribution to Z-boson decay flavor changing rates, Z qjqj', in- 
duced by radiative corrections from gluino-squark triangle diagrams of squarks flavor mixing, found p3t : 
Bjj' ~ IQ-^. This result is suspect since a more complete calculation of the effect performed subsequently 
obtained considerably smaller contributions. Both calculations rely on unrealistic inputs, including 
a wrong mass for the top-quark and too low values for the superpartners mass parameters. It is hoped 
that a complete updated study could be soon performed. In fact, during the last few years, the study 
of loop corrections in extended versions of the standard model has evolved into a streamlined activity. 
For instance, calculations of loop contributions to the magnetic moment of the r-lepton or of the heavy 
quarks, such as those reported in (two-Higgs doublets model) or in [ |65[ (minimal supersymmetric 
standard model) could be usefully transposed to the case of fermion pair production observables. 

The information from experimental searches on flavor changing physics at high energy colliders is 
rather meager | |6^ . Upper bounds for the leptonic Z-boson branching ratios, Bjji, are reported [ |67) at, 
1.7 10^^ for {efi + pe), 9.8 10~^ for (er 4- f e) and 1.7 10^^ for {fiT + f^i). No results have been quoted so 
far for d— or u— quark pairs production, reflecting the hard experimental problems faced in identifying 
quarks flavors at high energies. The prospect for experimental measurements at the future leptonic 
colliders is brightest for cases involving one top-quark owing to the easier kinematical identification 
offered by the large mass disparity in the final state jets. For leptonic colliders at energies above those of 
LEP, the reactions involving the production of Higgs or heavy Z'-gauge bosons which subsequently decay 
to fermion pairs could be effective sources of flavor non-diagonal effects, especially when a top-quark is 
produced. At still higher energies, in the TeV regime, the production subprocesses involving collisions of 
gauge bosons pairs radiated by the incident leptons, as in l~ -\- ^ W~ + W'^ + i' + D, could lead to 
flavor non-diagonal final states, such as, + P + t + c with rates of order a few fbarns |]68| . 



3.2 Choices of parameters and models 

Our main assumption in this work is that no other sources besides the R parity odd interactions contribute 
significantly to the flavor changing rates and CP asymmetries. However, to infer useful information 
from possible future experimental results, we must deal with two main types of uncertainties. The first 
concerns the family structure of the coupling constants. On this issue, one can only postulate specific 
hypotheses or make model-dependent statements. At this point, we may note that the experimental 
indirect upper bounds on single coupling constants are typically, A < 0.05 or A' < 0.05 times jq^^^, 
except for three special cases where strong bounds exist: A'^ < 3.9 10~"^( ioo'Gev )^( iooGey (0i^/3/3— 
decay H) A^g < 2 lO'^ mass H) and A^„, < 2. 10-2(^^^), /c = 1,2,3; m = 1,2], 
{K TTvv Q)- Strong bounds have been derived for products of couphng constants pairs in specific 
family configurations. For instance, a valuable source for the \ijk coupling constants is provided by the 
rare decays, + e~ + Cp |^9j, which probe the combinations of coupling constants, Fated = 

Ei(^7^7^)^^»fc^icrf- Except for the strong bound, Ff^^a + -Ffiii < 4.3 lO^^^^ [n 3e] the other 
combinations of coupling constants involving the third generation are less strongly bound, as for instance, 
-^1113 +-^3111 < 3.1 10^5 [t 3e] [^. Another useful source is provided by the neutrinoless double beta 
decay process ^ . The strongest bounds occur for the following configurations of flavour indices 

(using the reference value m = lOOGe^): A'ngA'jai < 7.9x10"^, X[i2K2i < 2.3x10"^ X[\^ <4.6xl0-^ 
quoting from p3] where the initial analysis of ||23| was updated. Finally, the strongest bounds deduced 
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from neutral mesons {BB, KK) mixing parameters are: F[^^i < 2 10 '"^i^iaai < 3.3 10 ^{221 < 
4.5 10-9, |1|, where F'^,^, = 1MJ^)2a^^^a:;^. 

The second type of uncertainties concerns the spectrum of scalar superpartners. At one extreme, are 
the experimental lower bounds, which reach for sleptons, 40 — 65 GeV, and for squarks, 90 — 200 GeV, 
and at the other extreme, the theoretical naturalness requirement which sets an upper bound at 1 TeV. 

In order to estimate the uncertainties on predictions emanating from the above two sources, it is 
necessary to delineate the dependence of amplitudes on sfermion masses. Examining the structure of 
the relevant contributions to flavor changing rates for, say, the lepton case, we note that the t-channel 
exchange tree amplitudes are given by a onefold summation over sfermions families, i^i" 
volving the combination of coupling constants, V'jj, = AiijA*]^j/. The typical structure for the leptonic 
loop amplitudes is a twofold summation over fermions and sfermions families, '^ij ljj>Fjj,{s + ie), 
where l^j j, = ^ijj^ijjr, and the loop integrals, F^j,, have a non-trivial dependence on the fermions and 



sfermions masses, as exhibited on the formulas derived in subsections 2.1, p.2| and 2.3 [see, e.g., eq.(|21|)]. 



The effective dependence on the superparticle masses involves ratios, mj/fh or s/fn , in such a way 
that the dependence is suppressed for large fh. (Obviously, s = m| for Z-boson pole observables.) In 
applications such as ours where, s > m^, all the fermions, with the exception of the top-quark, can be 
regarded as being massless. In particular, the first two light families (for either d, u) should have 
comparable contributions, the third family behaving most distinctly in the top-quark case. A quick 
analysis, taking the explicit mass factors into account, indicates that loop amplitudes should scale with 
sfermions masses as, (s/m^)", with a variable exponent ranging in the interval, 1 < n < 2. Any possible 
enhancement effect from the explicit sfermions mass factors in eq.(pl|) is moderated in the full result 
by the fact that the accompanying loop integral factor has itself a power decrease with increasing fh?'. 
Thus, the Z-boson pole rates should depend on the masses fh roughly as (1/m^)^", while the off Z-boson 
pole rates, being determined by the tree amplitudes, should behave more nearly as (l/m^)^. As for the 
asymmetries, since these are given by ratios of squared amplitudes, one expects them to have a weak 
sensitivity on the sfermion masses. 

To infer the physical implications on the RPV coupling constants, we avoid making too detailed model- 
dependent assumptions on the scalar superpartners spectrum. Thus, we shall neglect mass splittings 
between all the sfermions and set uniformly all sleptons, sneutrinos and squarks masses at a unique 
family (species) independent value, m, chosen to vary in the wide variation interval, 100 < to < 1000 
GeV. This prescription should suffice for the kind of semi-realistic predictions at which we are aiming. 
This approximation makes more transparent the dependence on the RPV coupling constants, which 
then involves the quadratic products designated by tjj, (tree) and l^jj, (loop), where the dummy family 
indices refer to sfermions (tree) and fermion-sfermions (loop). For off Z-boson-pole observables, flavor 
non diagonal rates are controlled by products of two different couplings, and asymmetries by 

normalized products of four different couplings, Im{tjy,rjj,)/\tjj,\'^. For Z-boson pole observables, rates 
and asymmetries are again controlled by products of two and four different coupling constants, l^'/j/p 
and Im{Vj''j*V'jj,)/\V'jj, p, respectively. Let us note that if the off-diagonal rates were dominated by some 
alternative mechanism, the asymmetries would then involve products of four different coupling constants 
rather than the above ratio. 

It is useful here to set up a catalog of the species and families configurations for the sfermions 
(tree) or fermion-sfermions (loop) involved in the various cases. In the tree level amplitudes, these 
configurations are for leptons: tjj, — i>iL (£-type), t^jj, = XiijX*ij,, ViL (7?.-type); for d- 

quarks, t-'j j, — X[jjX*jj,, ujl ("T^-type); for u-quarks, tjj, = X^jf.\\j,j^, dkR (£-type). In the loop level 
amplitudes, the coupling constants and internal fermion-sfermion configurations are for leptons: 

1% = [(4), (£;)]; Ifj' = Kj,>^U'k. (^) (r-type); 

^Jr^K,jK,j'^ [(^;)] (7^-type); 

for d-quarks: 

1%'->:u,Kj'u, [(^), (£)] ('C-type); 

^]r = >^^rK,j. lUl), Ul); (4;)] (7^-type); 

for u-quarks, 

vfr-y^uuKj'u^ [Gt), (£)] (/:-type). 

We shall present numerical results for a subset of the above list of cases. For leptons and d-quarks. 
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we shall restrict consideration to the 7?,- type terms which contribute to the Z-boson vertex function, Ar. 
We also retain the sleptons-quarks internal states for d-quark production (involving \ and the 

sleptons-leptons for lepton production (involving For the up-quark production, we consider the 

jC-type terms (involving and, for the off Z-boson pole case, omit the term ipLR in eq.(23) in 

the numerical results. 

Since the running family index in the parameters relevant to tree level amplitudes refers to sfermions, 
consistently with the approximation of a uniform family independent mass spectrum, we may as well 
consider that index as being fixed. Accordingly, we shall set these parameters at the reference value, 
tjj, = 10^^. In contrast to the off Z-boson pole rates, the asymmetries depend non trivially on the 
fermion mass spectrum through one of the two family indices in Vj j, (i or j) associated to fermions. 
To discuss our predictions, rather than going through the list of four distinct coupling constants, we 
shall make certain general hypotheses regarding the generation dependence of the RPV interactions for 
the fermionic index. At one extreme is the case where all three generations are treated alike, the other 
extreme being the case where only one generation dominates. We shall consider three different cases 
which are distinguished by the interval over which the fermions family indices are allowed to range in the 
quantities, Ijj,- We define Case A by the prescription of equal values for all three families of fermions 
(i — 1,2,3); Case B, for the second and third families (i = 2,3); and Case C, for the third family only 
(i = 3). For all three cases, we set the relevant parameters uniformly at the reference values, tjj, = 10^^. 
While the results in Case C reflect directly on the situation associated with the hypothesis of dominant 
third family configurations, the corresponding results in situations where the first or second family are 
assumed dominant, can be deduced by taking the differences between the results in Cases A and B and 
Cases B and C, respectively. 

In order to obtain non- vanishing CP asymmetries, we still need to specify a prescription for introducing 
a relative CP-odd complex phase, denoted ip, between the various RPV coupling constants. We shall 
set this at the reference value, ip = 7r/2. Since the CP asymmetries are proportional to the imaginary 
part of the phase factor, the requisite dependence is simply reinstated by inserting the overall factor, 
sin-f/). Different prescriptions must be implemented depending on whether one considers observables 
at or off the Z-boson pole. The Z-boson pole asymmetries are controlled by a relative complex phase 
between the combinations of coupling constants denoted, Zjj, only. For definiteness, we choose here to 
assign a non- vanishing complex phase only to the third fermion family, namely, arg{iy j,) = [0, 0, 7r/2], for 
[i or j — 1, 2, 3]. In fact, a relative phase between light families only would contribute insignificantly to the 
Z-boson pole asymmetries, because of the antisymmetry in a — > a' in eq. ([lo|) and the fact that F^'^ {^%) 
are approximately equal when the fermion index in a = (j,i) belongs to the two first families. The off 
Z-boson pole asymmetries are controlled by a relative complex phase between the tree and loop level 
amplitudes. For definiteness, we choose here to assign a vanishing argument to the coupling constants 
combination, t^j, appearing at tree level and non-vanishing arguments to the full set of loop amplitude 
combinations, namely, argltjj,) = 7r/2, where the fermion index {i or j as the case may be) varies over 
the ranges relevant to each of the three cases A, B, C. 



3.3 Numerical results and discussion 
3.3.1 Z-boson decays observables 

We start by presenting the numerical results for the integrated rates associated with Z-boson decays into 
fermion pairs. These are given for the d-quarks, leptons and u-quarks cases in Table |l]. We observe a 
fast decrease of rates with increasing values of the mass parameter, m. Our results can be approximately 
fitted by a power law dependence which is intermediate between m^^ and rfi^^. Explicitly, the Z-boson 
flavor non diagonal decay rates to d-quarks, leptons and u-quarks, are found to scale approximately as, 
Bjj' « ( ^"001''' )^( ^°°^''^ )^'^ X 10"^[5., 1., 2.], respectively. When a top-quark intermediate state is 
allowed in the loop amplitude, this dominates over the contributions from the light families. This is 
clearly seen on the d-quarks results which are somewhat larger than those for up-quarks and significantly 
larger than those for leptons, the more so for larger rh. This result is explained partly by the color factor, 
partly by the presence of the top-quark contribution only for the down-quarks case. For contributions 
involving other intermediate states than up-quarks, whether the internal fermion generation index in the 
RPV coupling constants, Xijk, runs over all three generations (Case A), the second and third generations 
(Case B) or the third generation only (Case C), we find that rates get reduced by factors roughly less than 
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2 in each of these stages. Therefore, this comparison indicates a certain degree of family independence 
for the Z-boson branching fractions for the cases where either leptons or d-quarks propagate inside the 
loops. 

Proceeding next to the CP-odd asymmetries, since these are proportional to ratios of the RPV coupling 
constants, it follows in our prescription of using uniform values for these, that asymmetries must be 
independent of the specific reference value chosen. As for their dependence on to, we see on Table |l] 
that this is rather strong and that the sense of variation with increasing ifi corresponds (for absolute 
values of Ajji) to a decrease for d-quarks and an increase for u-quarks and leptons. The comparison 
of different production cases shows that the CP asymmetries are largest, 0(10""'^), for d-quarks at small 
TO ~ lOOGey, and for u-quarks at large to ~ lOOOGey. For leptons, the asymmetries are systematically 
small, O(10~'^ — 10^^). The above features are explained by the occurrence for d-quarks production of 
an intermediate top-quark contribution and also by the larger values of the rates at large to in this case. 
The comparison of results in Cases A and B indicates that the first two light families give roughly equal 
contributions in all cases. 

For Case C, the CP-odd asymmetries are vanishingly small, as expected from our prescription of 
assigning the CP-odd phase, since Case C corresponds then to a situation where only single pairs of 
coupling constants dominate. Recall that for the specific cases considered in the numerical applications, 
namely, 7?.- type for d-quarks and leptons and >C-type for u-quarks, the relevant products of RPV cou- 
pling constants are, Kjj'^^ijJ^ \*jkKj'ki respectively, where the fermions generation index 
amongst the dummy indices pairs, (zj), (ik), refers to the third family. Non vanishing contributions to 
Ajj' could arise in Case C if one assumed that two pairs of the above coupling constants products with 
different sfcrmions indices dominate, and further requiring that these sfermions are not mass degener- 
ate. Another interesting possibility is by assuming that the hypothesis of single pair of RPV coupling 
constants dominance applies for the fields current basis. Applying then to the quark superfields the 
transformation matrices relating these to mass basis fields, say, in the distinguished choice where 
the flavor changing effects bear on u-quarks, amounts to perform the substitution, A-^j, — + \nkKlj^ 
where V is the CKM matrix. The CP-odd factor, for the d-quark case, say, acquires then the form, 
^ |A:f^,A:^t,|2/TO((yt)„^.(yt);;^.(yt);;^.,(yt)„^^.,), where the second factor on the right- 
hand side is recognized as the familiar plaquette term, proportional to the products of sines of all the 
CKM rotation angles times that of the CP-odd phase. 

It may be useful to examine the bounds on the RPV coupling constants implied by the current 
experimental limits on the flavor non diagonal leptonic widths |^, Bj^f, < [1.7, 9.8, 17.] 10~^ for the 
family couples, [J J' = 12, 23, 13]. The contributions associated with the A interactions can be directly 
deduced from the results in Table 0. Choosing the value, rh = 100 GeV, and writing our numerical result 

A ■ - J A* 

as, Bjj' w ( '^0 01^^' )^'^ 10~^, then under the hypothesis of a pair of dominant coupling constants, one 
deduces, Ai^jA* j, < [0.46, 1.1, 1.4], for all fixed choices of the family couples, i, j. (An extra factor 2 
in Bjji has been included to account for the antisymmetry property of \ijk-) For the A' interactions, 
stronger bounds obtain because of the extra color factor and of the internal top-quark contributions. A 
numerical calculation (not reported in Table [^) performed with the choice, to = lOOGeV^ for Case C, 
a'* a' 

gives us: Bjji « {^^^^-^^Yl.YI 10^^, which, by comparison with the experimental limits, yields the 
bounds: ^jjf^^'j,jk < [0.38, 0.91, 1.2]10^^, for the same family configurations, [J J' = 12, 23, 13], as 
above. These results agree in size to within a factor of 2 with results reported in a recently published 
work 

3.3.2 Fermion anti-fermion pair production rates 

Let us now proceed to the off Z-boson pole observables. The numerical results for the flavor non-diagonal 
integrated cross sections and CP asymmetries are shown in Table |l| for two selected values of the center of 
mass energy, y/s — 200 and 500 GeV. The numerical results displaying the variation of these observables 
with the center of mass of energy (fixed to) and with the superpartners mass parameter (fixed ^/s) are 
given in Fig.^ and Fig.^, respectively. All the results presented in this work include both photon and 
Z-boson exchange contributions. We observe here that the predictions for asymmetries are sensitive to 
the interference effects between photon and Z-boson exchange contributions. 

We discuss first the predictions for flavor non diagonal rates. We observe a strong decrease with 
increasing values of rh and a slow decrease with increasing values of ^/s. Following a rapid initial rise 
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Table 1: Flavor changing rates and CP asymmetries for d-quarks, leptons and u-quarks pair production 

in the three cases, appearing in hne entries as Cases A, B and C, which correspond to internal hnes 
belonging to all three families, the second and third families and the third family, respectively. The 
results for d-quarks and leptons, unlike those for up-quarks, are obtained in the approximation where 
one neglects the final fermions masses. The first four column fields (Z-pole column entry) show results 
for the Z-boson pole branching fractions Bjjr and asymmetries, Ajj'. The last four column fields (off 
Z-pole column entry) show results for the flavor non-diagonal cross sections, ajj', in fbarns and for the 
asymmetries, Ajj', with photon and Z-boson exchanges added in. The results in the two lines for the off 
Z-boson pole are associated to the two values for the center of mass energy, s^/-^ = 200, 500 GeV. The 
columns subentries indicated by m correspond to the sfermions mass parameter, m = 100, 1000 GeV. 
The notation d — n stands for 10"". 
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at threshold, the rates settle at values ranging between (10 — 10^^) fbarns for a wide interval of m 
values. The dependence on rh can be approximately represented as, ^ '^J ~ (1 — 10) 

fbarns ~ )^(10^^ — 1). The rate of decrease of ctjj' with to slows down with increasing s. It is 

interesting to note that if we had considered here constant values of the product A (m/lOOGeF), rather 
than constant values of A, the power dependence of rates on rh would be such as to lead to interestingly 
enhanced rates at large rh. 

The marked differences exhibited by the results for lepton pair production, apparent on windows (c) 
and (d) in Figures I and I are due to our deliberate choice of adding the s-channel v pole term for the 
lepton case while omitting it for the d-quark case. The larger rates found for leptons as compared to d- 
quarks, in spite of the extra color factor present for d-quarks (recall that the fjfj' reactions rates 

for down-quarks and up-quarks pick up an extra color factor Nc with respect to those for leptons) is thus 
explained by the strong enhancement induced by adding in the sneutrino exchange contribution. This 
choice was made here for illustrative purposes, setting for orientation the relevant coupling constant at the 
value, Xijj' — 0.1. The z/ propagator pole was smoothed out by employing the familiar shifted propagator 
prescription, (s — m? + zTOpF^;)^^, while describing approximately the sneutrinos decay width in terms 

of the RPV contributions alone, namely, T{Di ^ + ij) — '{^^ ' and r(t>i + dj) = Nc — ig^- 

Proceeding next to the CP-odd asymmetries, we note that since these scale as a function of the 
RPV coupling constants as, Im{ljj,ljj*)/\tjj,\'^, our present predictions are independent of the uniform 
reference value assigned to these coupling constants. If the generational dependence of the RPV cou- 
pling constants were to exhibit strong hierarchies, this peculiar rational dependence could induce strong 
suppression or enhancement factors. 

The cusps in the dependence of Ajj' on y/s (Fig. ||) occur at values of the center of mass energy where 
one crosses thresholds for fermion-antifermion (for the energies under consideration, ti) pair production, 
y/s = 2m f, and scalar superpartners pair production, ^/s = 2rh. These are the thresholds for the 
processes, l~ + ff or + /+ f f *, at which the associated loop amplitudes acquire finite 

imaginary parts. Correspondingly, in the dependence of Ajj' on rh (Fig. ^ the cusps appear at rh — 
\/s/2. We note on the results that the ti contributions act to suppress the asymmetries whereas the 
//* contributions rather act to enhance them. Sufficiently beyond these two-particle thresholds, the 
asymmetries vary weakly with rh. A more rapid variation as a function of energy occurs in the leptons 
production case due to the addition there of the sneutrino pole contribution. 

The comparison of results for asymmetries in Cases A, B, C reflects on the dependence of loop 
integrals with respect to the internal fermions masses. An examination of Table |^ reveals that for leptons 
and up-quarks, where intermediate states involve leptons or d-quarks, all three families have nearly equal 
contributions. The results for down-quarks production are enhanced because of the intermediate top- 
quark contribution, which dominates over that of lighter families. However, this effect is depleted when 
the finite imaginary part from tt sets in. The asymmetries for up-quarks production assume values in 
the range, 10~^ — 10"'^, irrespective of the fact that the final fermions belong to light or heavy families. 

4 Conclusions 

The two-body production at high energy leptonic colliders of fermion pairs of different families could 
provide valuable information on the flavor structure of the R parity odd Yukawa interactions. One 
can only wish that an experimental identification of lepton and quark flavors at high energies becomes 
accessible in the future. Although the supersymmetric loop corrections to these processes may not 
be as strongly suppressed as their standard model counterparts, one expects that the degeneracy or 
alignment constraints on the scalar superpartners masses and flavor mixing should severely bound their 
contributions. Systematic studies of the supersymmetry corrections to the flavor changing rates and CP 
asymmetries in fermion pair production should be strongly encouraged. 

An important characteristic of the R parity odd interactions is that they can contribute to integrated 
rates at tree level and to CP asymmetries through interference terms between the tree and loop ampli- 
tudes. While we have restricted ourselves to the subset of loop contributions associated with Z-boson 
exchange, a large number of contributions, involving quark-sleptons or lepton-squarks intermediate states 
in various families configurations, could still occur. The contributions to rates and asymmetries depend 
strongly on the values of the R parity odd coupling constants. Only the rates are directly sensitive to the 
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supersymmetry breaking scale. To circumvent the uncertainties from the sparticles spectrum, we have 
resorted to the simphfying assumption that the scalar superpartners mass differences and mixings can 
be neglected. We have set the RPV coupling constants at a uniform value while sampling a set of cases 
from which one might reconstruct the family dependence of the RPV coupling constants. We have also 
embedded a CP complex phase in the RPV coupling constants in a specific way, meant to serve mainly 
as an illustrative example. Although the representative cases that we have considered represent a small 
fraction of the host of possible variations, they give a fair idea of the sizes to expect. Since these processes 
cover a wide range of family configurations, one optimitistic possibility could be that one specific entry 
for the family configurations would enter with a sizeable RPV coupling constant. 

The contributions to the fiavor changing rates have a strong sensitivity on the RPV coupling constants 
and the superpartners mass, involving high powers of these parameters. We find a generic dependence for 
the flavor changing Z-boson decay branching ratios of form, {■^^)'^{^)'^'^ 10~^. For the typical bounds 
on the RPV coupling constants, it appears that these branchings are three order of magnitudes below the 
current experimental sensitivity. At higher energies, the flavor changing rates are in order of magnitude, 
(er^)^('^)^ ~ (1 ~ 10) fbarns. Given the size for the typical integrated luminosity, C = 50/6~^/year, 
anticipated at the future leptonic machines, one can be moderately optimistic on the observation of clear 
signals. 

The Z-boson pole CP-odd asymmetries are of order, (10~^ — 10^^) sinip. For the off Z-boson pole 
reactions, a CP-odd phase, ip^ embedded in the RPV coupling constants shows up in asymmetries with 
reduced strength, (10^^ — 10^"^) sin?/' for leptons, d-quarks and u-quarks. The largely unknown structure 
of the RPV coupling constants in flavor space leaves room for good or bad surprises, since the peculiar 
rational dependence on the couphng constants, /m(AA*AA'^)/A'', and similarly with A A', may lead to 
strong enhancement or suppression factors. 
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Figure 3: Integrated flavor non-diagonal cross sections (left hand side windows) and asymmetries (right 
hand side windows) as functions of the center of mass energy in the production of down-quark-antiquark 
pairs (two upper figures (a) and (6)), lepton-antilepton pairs (two intermediate figures (c)and (d)) and 
up-quark-antiquark pairs of typo tc + ct (two lower figures (e) and (f)). The tree level amplitude includes 
only the t-channel contribution for the d-quark case, both t- and s-channel exchange contributions for the 
lepton case, and the u-channel exchange for the up-quark case. The one-loop amplitudes, with both Z- 
boson and photon exchange contributions, include all three internal fermions generations, corresponding 
to Case A. Three choices for the superpartners uniform mass parameter, m, are considered: lOOGeV 
(continuous lines), 200GeV (dashed-dotted lines), 500GeV (dashed lines). 
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Figure 4: Integrated flavor non-diagonal cross sections (left hand side windows) and CP-odd asymmetries 
(right hand side windows) as functions of the scalar superpartners mass parameter, m, in the production 
of down-quark-antiquark pairs (two upper figures (^f and (6)), lepton-antilepton pairs (two intermediate 

figures (c)and (d)) and up-quark-antiquark pairs of type tc or ct (two lower figures (e) and (f)). The 



